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Abstract 

This paper is dedicated to Oleg Viro on his 60-th birthday. In the paper we explore how the 
Potts model in statistical mechanics is related to Khovanov homology. This exploration is made 
possible because the underlying combinatorics for the bracket state sum for the Jones polynomial 
is shared by the Potts model for planar graphs. We show that Euler characteristics of Khovanov 
homology figure in the computation of the Potts model at certain imaginary temperatures and 
that these aspects of the Potts model can be reformulated as physical quantum amplitudes via 
Wick rotation. The paper concludes with a new conceptually transparent quantum algorithm for 
the Jones polynomial and with many further questions about Khovanov homology. 

1 Introduction 

This paper is about Khovanov homology and its relationships with finite combinatorial statistical 
mechanics models such as the Ising model and the Potts model. 

Partition functions in statistical mechanics take the form 
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where a runs over the different physical states of a system G and E{a) is the energy of the state 
a. The probability for the system to be in the state a is taken to be 
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Since Onsager's work, showing that the partition function of the Ising model has a phase transi- 
tion, it has been a significant subject in mathematical physics to study the properties of partition 



functions for simply defined models based on graphs G. The underlying physical system is mod- 
eled by a graph G and the states a are certain discrete labellings of G. The reader can consult 
Baxter's book 01 for many beautiful examples. The Potts model, discussed below, is a general- 
ization of the Ising model and it is an example of a statistical mechnanics model that is intimately 
related to knot theory and to the Jones polynomial J8J. Since there are many connections between 
statistical mechanics and the Jones polynomial, it is remarkable that there have not been many 
connections of statistical mechanics with categorifications of the Jones polynomial. The reader 
will be find other points of view in [[51 [6J. 

The partition function for the Potts model is given by the formula 

P G (Q,T) = £ e WwW =Y^e KE ^ 

where a is an assignment of one element of the set {1, 2, • • • , Q} to each node of the graph G, and 
E{a) denotes the number of edges of a graph G whose end-nodes receive the same assignment 
from a. In this model, a is regarded as a physical state of the Potts system and E(a) is the energy 
of this state. Here K = J-g^ where J is plus or minus one (ferromagnetic and anti-ferromagnetic 
cases), k is Boltzmann's constant and T is the temperature. The Potts partition function can be 
expressed in terms of the dichromatic polynomial of the graph G. Letting 



it is shown in lfT2l [T3l that the dichromatic polynomial Z[G](v, Q) for a plane graph can be 
expressed in terms of a bracket state summation of the form 

{X} = {X} + Q'HX} 

with 

{0} = Q*- 

Then the Potts partition function is given by the formula 

P G = Q N ' 2 {K{G)} 

where K{G) is an alternating link diagram associated with the plane graph G so that the projec- 
tion of K(G) to the plane is a medial diagram for the graph. This translation of the Potts model 
in terms of a bracket expansion makes it possible to examine how the Khovanov homology of the 
states of the bracket is related to the evaluation of the partition function. 

There are five sections in this paper beyond the introduction. In Section 2 we review the 
definition of Khovanov homology and observe, in parallel with G71 , that it is very natural to 
begin by defining the Khovanov chain complex via enhanced states of the bracket polynomial 
model for the Jones polynomial. In fact, we begin with Khovanov's rewrite of the bracket state 
sum in the form 

(X> = <X>-?00 
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with (O) = (q + We rewrite the state sum formula for this version of the bracket in terms 
of enhanced states (each loop is labeled +1 or —1 corresponding to q +1 and q^ 1 respectively) 
and show how, collecting terms, the formula for state sum has the form of a graded sum of Euler 
characteristics 

3 i 3 3 

where C u is a module with basis the set of enhanced states s with i smoothings of type B (see 
Section 2 for definitions), TP 3 is its homology and j = j(s) where 

j(s) = n B (s) + X(s) 

where tib(s) denotes the number of 5-smoothings in s and A(s) denotes the number of loops 
with positive label minus the number of loops with negative label in s. This formula suggests that 
there should be differentials d : C lJ — > C+ij so ^at j is preserved under the differential. We 
show that the restriction j(s) = j(d(s)) uniquely determines the differential in the complex, and 
how this leads to the Frobenius algebra structure that Khovanov used to define the differential. 
This part of our remarks is well-known, but I believe that the method by which we arrive at the 
graded Euler characteristic is particularly useful for our subsequent discussion. We see clearly 
here that one should look for subcomplexes on which Euler characteristics can be defined, and 
one should attempt to shape the state summation so that these characteristics appear in the state 
sum. This happens miraculously for the bracket state sum and makes the combinatorics of that 
model dovetail with the Khovanov homology of its states, so that the bracket polynomial is seen 
as a graded Euler characteric of the homology theory. The section ends with a discussion of how 
Grassmann algebra can be used, in analogy with de Rahm cohomology, to define the integral 
Khovanov chain complex. We further note that this analogy with de Rahm cohomology leads to 
other possibilities for chain complexes associated with the bracket states. These complexes will 
be the subject of a separate paper. 

In Section 3 we recall the definition of the dichromatic polynomial and the fact lfT2l that 
the dichromatic polynomial for a planar graph G can be expressed as a special bracket state 
summation on an associated alternating knot K(G). We recall that, for special values of its two 
parameters, the Potts model can be expressed as a dichromatic polynomial. This sets the stage for 
seeing what is the role of Khovanov homology in the evaluations of the Potts model and of the 
dichromatic polynomial. In both cases (dichromatic polynomial in general and the Potts model 
in particular) one finds that the state summation does not so easily rearrange itself as a sum over 
Euler characteristics, as we have explained in Section 2. The states in the bracket summation 
for the dichromatic polynomial of K(G) are the same as the states for the bracket polynomial of 
K(G), so the Khovanov chain complex is present at the level of the states. 

We clarify this relationship by using a two-variable bracket expansion, the p-bracket, that 
reduces to the Khovanov version of the bracket as a function of q when p is equal to one.S 

IX] = IX] - qplX] 
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with 

[0] = 3 + ? _1 - 

We can regard this expansion as an intermediary between the Potts model (dichromatic polyno- 
mial) and the topological bracket. The p-bracket can be rewritten in the following form: 

\K] = 52(-pyq j dim(C ij ) 

3 i j 

where we define C u to be the linear span of the set of enhanced states with rag(s) = i and 
j(s) = j. Then the number of such states is the dimension dim(C^). We have expressed this 
p-bracket expansion in terms of generalized Euler characteristics of the complexes C'-> : 

X P (C' j ) = -£(- P ydim(C«). 

i 

These generalized Euler characteristics become classical Euler characteristics when p — 1, and, 
in that case, are the same as the Euler characteristic of the homology. With p not equal to 1, we 
do not have direct access to the homology. In that case the polynomial is expressed in terms of 
ranks of chain modules, but not in terms of ranks of corresponding homology groups. 

In Section 3 we analyze those cases of the Potts model where p = 1 (so that Euler character- 
istics of Khovanov homology appear as coefficients in the Potts partition function) and we find 
that at criticality this requires Q = 4 and e K = — 1 (K = J/kT as in the second paragraph of 
this introduction), hence an imaginary value of the temperature. When we simply require that 
p — 1, not necessarily at criticality, then we find that for Q = 2 we have e K = ±1 For Q = 3, we 
have e K = ~ 1± ^ t , For Q = 4 we have e K = —1. For Q > 4 it is easy to verify that e K is real 
and negative. Thus in all cases of p = 1 we find that the Potts model has complex temperature 
values. Further work is called for to see how the evaluations of the Potts model at complex values 
influence its behaviour for real temperature values in relation to the Khovanov homology. Such 
relationships between complex and real temeperature evaluations already known for the accum- 
lation of the zeroes of the partition function (the Lee- Yang zeroes ll20l ). In section 5 we return to 
the matter of imaginary temperature. This section is described below. 

In Section 4 we discuss Stosic's categorification of the dichromatic polynomial. Stosic's cat- 
egorification involves using a differential motivated directly by the graphical structure, and gives 
a homology theory distinct from Khovanov homology. We examine Stosic's categorification in 
relation to the Potts model and again show that it will work (in the sense that the coefficients of of 
partition function are Euler characteristics of the homology) when the temperature is imaginary. 
Specifically, we find this behaviour for K = in + ln(q + q 2 + • • • + q n ) , and again the challenge is 
to find out what is the influence of this graph homology on the Potts model at real temperatures. 
The results of this section do not require planarity of the graph G. More generally, we see that in 
the case of this model, the differential in the homology is related to the combinatorial structure 
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of the physical model. In that model a given state has graphical regions of constant spin (calling 
the discrete assignments {1, 2, • • • , Q} the spins). We interpret the partial differentials in this 
model as taking two such regions and making them into a single region, by adding an edge that 
joins them, and re-assigning a spin to the new region that is a combination of the spins of the 
two formerly separate regions. This form of partial differential is a way to think about relating 
the different states of the physical system. It is however, not directly related to classical physical 
process since it invokes a global change in the spin confinguration of the model. (One might con- 
sider such global transitions in quantum processes.) The Stosic homology measures such global 
changes, and it is probably the non-classical physicality of this patterning that makes it manifest 
at imaginary temperatures. The direct physical transitions from state to state that are relevant 
to the classical physics of the model may indeed involve changes of the form of the regions of 
constant spin, but this will happen by a local change so that two disjoint regions with the same 
spin join into a single region, or a single region of constant spin bifurcates into two regions with 
this spin. These are the sort of local-classical-time transitions that one works with in a statisti- 
cal mechanics model. Such transitions are part of the larger and more global transitions that are 
described by the differentials in our interpretation of the Stosic homology for the Potts model. 
Obviously, much more work needs to be done in this field. We have made first steps in this paper. 

In Section 5 we formulate a version of Wick rotation for the Potts model so that it is seen 
(for imaginary temperature) as a quantum amplitude. In this way, for those cases where the 
temperature is pure imaginary, we obtain a quantum physical interpretation of the Potts model 
and hence a relationship of Khovanov homology with quantum amplitudes at the special values 
discussed in Section 3. In Section 6, we remark that the bracket state sum itself can be given a 
quantum statistical interpretation, by choosing a Hilbert space whose basis is the set of enhanced 
states of a diagram K . We use the evaluation of the bracket at each enhanced state as a matrix 
element for a linear transformation on the Hilbert space. This transformation is unitary when 
the bracket variable (here denoted as above by q) is on the unit circle. In this way, by using the 
Hadamard test, we obtain a new quantum algorithm for the Jones polynomial at all values of 
the Laurent polynomial variable that lie on the unit circle in the complex plane. This is not an 
efficient quantum algorithm, but it is conceptually transparent and it will allow us in subsequent 
work to analyze relationships between Khovanov homology and quantum computation. 

Acknowledgement. It gives the author of this paper great pleasure to acknowledge a helpful 
conversation with John Baez. 



2 Khovanov Homology 

In this section, we describe Khovanov homology along the lines of |fT9l [2l, and we tell the story 
so that the gradings and the structure of the differential emerge in a natural way. This approach to 
motivating the Khovanov homology uses elements of Khovanov's original approach, Viro's use of 
enhanced states for the bracket polynomial ll2~7ll . and Bar-Natan's emphasis on tangle cobordisms 
0. We use similar considerations in our paper ifTTll . 
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Two key motivating ideas are involved in finding the Khovanov invariant. First of all, one 
would like to categorify a link polynomial such as (K). There are many meanings to the term 
categorify, but here the quest is to find a way to express the link polynomial as a graded Euler 
characteristic (K) = XqCH{K)) for some homology theory associated with (K). 

The bracket polynomial [1 Q model for the Jones polynomial [|8l|9l[T0l|28l * s usually described 
by the expansion 

and we have 

(KO) = (-A 2 -A- 2 )(K) 
(V> = (-A- 3 )(w> 

Letting c(K) denote the number of crossings in the diagram K, if we replace (K) by A~ c ^ (K) , 
and then replace A 2 by — q~ l , the bracket will be rewritten in the following form: 

<X> = (X> -<?<)<> 

with (O) = (<? + It is useful to use this form of the bracket state sum for the sake of the 
grading in the Khovanov homology (to be described below). We shall continue to refer to the 
smoothings labeled q (or A^ 1 in the original bracket formulation) as B-smoothings. We should 
further note that we use the well-known convention of enhanced states where an enhanced state 
has a label of 1 or X on each of its component loops. We then regard the value of the loop q + q~ l 
as the sum of the value of a circle labeled with a 1 (the value is q) added to the value of a circle 
labeled with an X (the value is q~ l ). We could have chosen the more neutral labels of +1 and 
— 1 so that 

q +1 ^+1^1 

and 

q' 1 -1 X, 

but, since an algebra involving 1 and X naturally appears later, we take this form of labeling from 
the beginning. 

To see how the Khovanov grading arises, consider the form of the expansion of this version of 
the bracket polynonmial in enhanced states. We have the formula as a sum over enhanced states 

s : 

(K) = J2{-l) nB{s) q j{s) 

s 

where tib(s) is the number of 5-type smoothings in s, A(s) is the number of loops in s labeled 
1 minus the number of loops labeled X, and j(s) = n B (s) + A(s). This can be rewritten in the 
following form: 

(A') = 53(-l)Y'dim(C ij ') 

* i5 
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where we define C y to be the linear span (over k = Z/2Z as we will work with mod 2 coeffi- 
cients) of the set of enhanced states with %(s) = i and j(s) = j. Then the number of such states 
is the dimension dim{C^). 

We would like to have a bigraded complex composed of the C lJ with a differential 

d . c 3 — >C i+lj . 

The differential should increase the homological grading i by 1 and preserve the quantum grading 
j. Then we could write 

( K ) = £^£(-i)^™(n =J2 ( i j x(c j ), 

j i j 

where ) is the Euler characteristic of the subcomplex C*- 7 for a fixed value of j. 

This formula would constitute a categorification of the bracket polynomial. Below, we shall see 
how the original Khovanov differential d is uniquely determined by the restriction that j(ds) = 
j(s)for each enhanced state s. Since j is preserved by the differential, these subcomplexes C' j 
have their own Euler characteristics and homology. We have 

X (H(C' j )) = x(C' 3 ) 
where HiC' '-*) denotes the homology of the complex C*K We can write 

3 

The last formula expresses the bracket polynomial as a graded Euler characteristic of a homology 
theory associated with the enhanced states of the bracket state summation. This is the categorifi- 
cation of the bracket polynomial. Khovanov proves that this homology theory is an invariant of 
knots and links (via the Reidemeister moves of Figure 1), creating a new and stronger invariant 
than the original Jones polynomial. 

We will construct the differential in this complex first for mod-2 coefficients. The differential 
is based on regarding two states as adjacent if one differs from the other by a single smoothing 
at some site. Thus if (s, r) denotes a pair consisting in an enhanced state s and site r of that 
state with r of type A, then we consider all enhanced states s' obtained from s by smoothing at 
r and relabeling only those loops that are affected by the resmoothing. Call this set of enhanced 
states S'[s, r]. Then we shall define the partial differential d T (s) as a sum over certain elements 
in S'[s, r], and the differential by the formula 

a( a ) = 5X*) 

T 

with the sum over all type A sites r in s. It then remains to see what are the possibilities for d r (s) 
so that j(s) is preserved. 
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Figure 1 : Reidemeister Moves 



Note that if s' E S'[ 



s, t), then wb(s') = Ub(s) + 1. Thus 



jV) = n B {s') + A(s') = 1 + n B (s) + A(s')- 



From this we conclude that j(s) = j(s 



') if and only if A(s') 



A(s) - 1. Recall that 



\(s) = [s : +] - [s : 



where [s : +] is the number of loops in s labeled +1, [s : — ] is the number of loops labeled — 1 
(same as labeled with X) and j(s) = n B (s) + A(s). 

Proposition. The partial differentials d T (s) are uniquely determined by the condition that j(s') = 
j(s) for all s' involved in the action of the partial differential on the enhanced state s. This unique 
form of the partial differential can be described by the following structures of multiplication and 
comultiplication on the algebra A = k[X]/(X 2 ) where k = Z/2Z for mod-2 coefficients, or 
k = Z for integral coefficients. 

1. The element 1 is a multiplicative unit and X 2 = 0. 

2. A(l) = 1 <g> X + X ® 1 and A(X) = X®X. 

These rules describe the local relabeling process for loops in a state. Multiplication corresponds 
to the case where two loops merge to a single loop, while comultiplication corresponds to the 
case where one loop bifurcates into two loops. 

Proof. Using the above description of the differential, suppose that there are two loops at r that 
merge in the smoothing. If both loops are labeled 1 in s then the local contribution to A(s) is 2. 
Let s' denote a smoothing in S[s, r]. In order for the local A contribution to become 1, we see 
that the merged loop must be labeled 1. Similarly if the two loops are labeled 1 and X, then the 
merged loop must be labeled X so that the local contribution for A goes from to —1. Finally, 
if the two loops are labeled X and X, then there is no label available for a single loop that will 
give —3, so we define d to be zero in this case. We can summarize the result by saying that there 
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is a multiplicative structure m such that m(l, 1) = 1, m(l, X) = m(X, 1) = x, m(X, X) = 0, 
and this multiplication describes the structure of the partial differential when two loops merge. 
Since this is the multiplicative structure of the algebra A = k[X]/ (X 2 ), we take this algebra as 
summarizing the differential. 

Now consider the case where s has a single loop at the site r. Smoothing produces two loops. 
If the single loop is labeled X, then we must label each of the two loops by X in order to make A 
decrease by 1. If the single loop is labeled 1, then we can label the two loops by X and 1 in either 
order. In this second case we take the partial differential of s to be the sum of these two labeled 
states. This structure can be described by taking a coproduct structure with A(X) = X <g> X and 
A(l) = 1®X+X®1. We now have the algebra ,4 = k[X]/(X 2 ) with product m : A® A — ► A 
and coproduct A : A — > A® A, describing the differential completely. This completes the 
proof. 4(t 

Partial differentials are defined on each enhanced state s and a site r of type A in that state. 
We consider states obtained from the given state by smoothing the given site r. The result of 
smoothing r is to produce a new state s' with one more site of type B than s. Forming s' from 
s we either amalgamate two loops to a single loop at r, or we divide a loop at r into two dis- 
tinct loops. In the case of amalgamation, the new state s acquires the label on the amalgamated 
circle that is the product of the labels on the two circles that are its ancestors in s. This case of 
the partial differential is described by the multiplication in the algebra. If one circle becomes 
two circles, then we apply the coproduct. Thus if the circle is labeled X, then the resultant two 
circles are each labeled X corresponding to A(X) = X <g> X. If the orginal circle is labeled 
1 then we take the partial boundary to be a sum of two enhanced states with labels 1 and X in 
one case, and labels X and 1 in the other case, on the respective circles. This corresponds to 
A(l) = 1 ® X + X ® 1. Modulo two, the boundary of an enhanced state is the sum, over all 
sites of type A in the state, of the partial boundaries at these sites. It is not hard to verify directly 
that the square of the boundary mapping is zero (this is the identity of mixed partials!) and that 
it behaves as advertised, keeping j(s) constant. There is more to say about the nature of this 
construction with respect to Frobenius algebras and tangle cobordisms. In Figures 2 and 3 we 
illustrate how the partial boundaries can be conceptualized in terms of surface cobordisms. The 
equality of mixed paritals corresponds to topological equivalence of the corresponding surface 
cobordisms, and to the relationships between Frobenius algebras and the surface cobordism cat- 
egory. The proof of invariance of Khovanov homology with respect to the Reidemeister moves 
(respecting grading changes) will not be given here. See [HH |2l |3). It is remarkable that this 
version of Khovanov homology is uniquely specified by natural ideas about adjacency of states 
in the bracket polynomial. 

Remark on Integral Differentials. Choose an ordering for the crossings in the link diagram K 
and denote them by 1, 2, • ■ ■ n. Let s be any enhanced state of K and let 9j(s) denote the chain 
obtained from s by applying a partial boundary at the i-th site of s. If the i-th site is a smoothing 
of type A' 1 , then di(s) = 0. If the i-th site is a smoothing of type A, then di(s) is given by the 
rules discussed above (with the same signs). The compatibility conditions that we have discussed 
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Figure 2: SaddlePoints and State Smoothings 




F G H 

Figure 3: Surface Cobordisms 
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show that partials commute in the sense that dt (dj(s)) = dj (di(s)) for all i and j. One then defines 
signed boundary formulas in the usual way of algebraic topology. One way to think of this regards 
the complex as the analogue of a complex in de Rahm cohomology. Let {dx\, dx 2 , ■ ■ • , dx n ] be 
a formal basis for a Grassmann algebra so that dxi A dxj = —dxi A dxj Starting with enhanced 
states s in C°(K) (that is, state with all A-type smoothings) Define formally, di(s) = di(s)dxi 
and regard di(s) as identical with di(s) as we have previously regarded it in C X {K). In general, 
given an enhanced state s in C k (K) with £>-smoothings at locations i\ < i 2 < ■ ■ ■ < ik, we 
represent this chain as s dx^ A ■ ■ ■ A dx ik and define 

n 

d(s dxi ± A ■ ■ ■ A dx ik ) = <9?(s) dxj A dx^ A ■ ■ ■ A dxi k , 

i=i 

just as in a de Rahm complex. The Grassmann algebra automatically computes the correct signs 
in the chain complex, and this boundary formula gives the original boundary formula when we 
take coefficients modulo two. Note, that in this formalism, partial differentials di of enhanced 
states with a £>-smoothing at the site i are zero due to the fact that dxi A dxi = in the Grass- 
mann algebra. There is more to discuss about the use of Grassmann algebra in this context. For 
example, this approach clarifies parts of the construction in lfT8l . 

It of interest to examine this analogy between the Khovanov (co)homology and de Rahm 
cohomology. In that analogy the enhanced states correspond to the differentiable functions on a 
manifold. The Khovanov complex C k (K) is generated by elements of the form s dx^ A • • • A dx ik 
where the enhanced state s has £>-smoothings at exactly the sites If we were to follow 

the analogy with de Rahm cohomology literally, we would define a new complex DR(K) where 
DR k (K) is generated by elements s dx i± A • ■ ■ A dx ik where s is any enhanced state of the link 
K. The partial boundaries are defined in the same way as before and the global boundary formula 
is just as we have written it above. This gives a new chain complex associated with the link K . 
Whether its homology contains new topological information about the link K will be the subject 
of a subsequent paper. 

A further remark on de Rham cohomology. There is another deep relation with the de Rham 
complex: In [j2TTl it was observed that Khovanov homology is related to Hochschild homology 
and Hochschild homology is thought to be an algebraic version of de Rham chain complex (cyclic 
cohomology corresponds to de Rham cohomology), compare [22]. 



3 The Dichromatic Polynomial and the Potts Model 

We define the dichromatic polynomial as follows: 

Z[G](v, Q) = Z[G'](v, Q) + vZ[G"](v, Q) 
Z[*UG] =QZ[G]. 
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where G' is the result of deleting an edge from G, while G" is the result of contracting that 
same edge so that its end-nodes have been collapsed to a single node. In the second equation, • 
represents a graph with one node and no edges, and • U G represents the disjoint union of the 
single-node graph with the graph G. 

In HHEl it is shown that the dichromatic polynomial Z[G](v, Q) for a plane graph can be 
expressed in terms of a bracket state summation of the form 

{X} = {X} + 0"'«OO 

with 

{Ctf = Q*. 

Here 

Z[G](v 1 Q)=Q N / 2 {K(G)} 

where K(G) is an alternating link diagram associated with the plane graph G so that the projec- 
tion of K(G) to the plane is a medial diagram for the graph. Here we use the opposite convention 
from [13 in associating crossings to edges in the graph. We set K(G) so that smoothing K(G) 
along edges of the graph give rise to £>-smoofhings of K(G) . See Figure 4. The formula above, in 
bracket expansion form, is derived from the graphical contraction-deletion formula by translating 
first to the medial graph as indicated in the formulas below: 

z[X} = z[X} + vzOQ. 

Z[RUK] = QZ[K}. 

Here the shaded medial graph is indicated by the shaded glyphs in these formulas. The medial 
graph is obtained by placing a crossing at each edge of G and then connecting all these crossings 
around each face of G as shown in Figure 4. The medial can be checkerboard shaded in relation 
to the original graph G (this is usually called the Tait Checkerboard Graph after Peter Guthrie Tait 
who introduced these ideas into graph theory), and encoded with a crossing structure so that it 
represents a link diagram. R denotes a connected shaded region in the shaded medial graph. Such 
a region corresponds to a collection of nodes in the original graph, all labeled with the same color. 
The proof of the formula Z\G] = Q N ^ 2 {K(G)} then involves recounting boundaries of regions 
in correspondence with the loops in the link diagram. The advantage of the bracket expansion 
of the dichromatic polynomial is that it shows that this graph invariant is part of a family of 
polynomials that includes the Jones polynomial and it shows how the dichromatic polynomial for 
a graph whose medial is a braid closure can be expressed in terms of the Temperley-Lieb algebra. 
This in turn reflects on the sturcture of the Potts model for planar graphs, as we remark below. 

It is well-known that the partition function Pg(Q, T) for the Q-state Potts model in statistical 
mechanics on a graph G is equal to the dichromatic polynomial when 

T 1 

v = e ^ —i 
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Figure 4: Medial Graph, Tait Checkerboard Graph and K(G) 



where T is the temperature for the model and k is Boltzmann's constant. Here J = ±1 accord- 
ing as we work with the ferromagnetic or anti-ferromagnetic models (see (U Chapter 12). For 
simplicity we denote 

* = 4; 

kT 



so that 



We have the identity 



v = e K - 1. 



P G (Q,T) = Z[G](e K -1,Q). 
The partition function is given by the formula 

Pg(Q,T) = £e^) 

where a is an assignment of one element of the set {1, 2, • • • , Q} to each node of the graph G, and 
E{a) denotes the number of edges of the graph G whose end-nodes receive the same assignment 
from a. In this model, a is regarded as a physical state of the Potts system and E{a) is the energy 
of this state. Thus we have a link diagrammatic formulation for the Potts partition function for 
planar graphs G. 

P G (Q,T) = Q N / 2 {K(G)}(Q lV = e K - 1) 
where N is the number of nodes in the graph G. 

This bracket expansion for the Potts model is very useful in thinking about the physical struc- 
ture of the model. For example, since the bracket expansion can be expressed in terms of the 
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Temperley-Lieb algebra one can use this formalism to express the expansion of the Potts model 
in terms of the Temperley-Lieb algebra. This method clarifies the fundamental relationship of 
the Potts model and the algebra of Temperley and Lieb. Furthermore the conjectured critical 
temperature for the Potts model occurs for T when Q~^v = 1. We see clearly in the bracket 
expansion that this value of T corresponds to a point of symmetry of the model where the value 
of the partition function does not depend upon the designation of over and undercros sings in the 
associated knot or link. This corresponds to a symmetry between the plane graph G and its dual. 

We first analyze how our heuristics leading to the Khovanov homology looks when general- 
ized to the context of the dichromatic polynomial. (This is a different approach to the question 
than the method of Stosic ll26l or ll24l [71, but see the next section for a discussion of Stosic's 
approach to categorifying the dichromatic polynomial.) We then ask questions about the rela- 
tionship of Khovanov homology and the Potts model. It is natural to ask such questions since the 
adjacency of states in the Khovanov homology corresponds to an adjacency for energetic states 
of the physical system described by the Potts model, as we shall describe below. 

For this purpose we now adopt yet another bracket expansion as indicated below. We call this 
two-variable bracket expansion the p-bracket. It reduces to the Khovanov version of the bracket 
as a function of q when p is equal to one. 

[X] = [X] - <?pD0 

with 

[0] = 5 + 5" 1 - 

We can regard this expansion as an intermediary between the Potts model (dichromatic polyno- 
mial) and the topological bracket. When p = 1 we have the topological bracket expansion in 
Khovanov form. When 

-qp = Cf^v 

and 

q + q' 1 = 

we have the Potts model. We shall return to these parametrizations shortly. 
Just as in the last section, we have 

[K] = J2(-p) nB{s) q j{s) 

s 

where n B (s) is the number of 5-type smoothings in s, A(s) is the number of loops in s labeled 
1 minus the number of loops labeled X, and j(s) = n B (s) + A(s). This can be rewritten in the 
following form: 

[K] = J2(- P yq j dim(C ij ) 

3 i j 
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where we define C'- J to be the linear span of the set of enhanced states with ub{s) = i and 
j(s) = j. Then the number of such states is the dimension dim{C^). Now we have expressed this 
general bracket expansion in terms of generalized Euler characteristics of the complexes: 

x P (cn = Y,(- P ydim(c ij ). 

i 

These generalized Euler characteristics become classical Euler characteristics when p — 1, and, 
in that case, are the same as the Euler characteristic of the homology. With p not equal to 1, we 
do not have direct access to the homology. 

Nevertheless, I believe that this raises a significant question about the relationship of [K] (q, p) 
with Khovanov homology. We get the Khovanov version of the bracket polynomial for p = 1 so 
that for p = 1 we have 

[K](q,p = 1) = J2q j Xp (C- j ) =E^(C ,J )). 

3 3 

Away from p — 1 one can ask what is the influence of the homology groups on the coefficients 
of the expansion of p), and the corresponding questions about the Potts model. This is a 

way to generalize questions about the relationship of the Jones polynomial with the Potts model. 
In the case of the Khovanov formalism, we have the same structure of the states and the same 
homology theory for the states in both cases, but in the case of the Jones polynomial (p-bracket 
expansion with p — 1) we have expressions for the coefficients of the Jones polynomial in terms 
of ranks of the Khovanov homology groups. Only the ranks of the chain complexes figure in the 
Potts model itself. Thus we are suggesting here that it is worth asking about the relationship of the 
Khovanov homology with the dichromatic polynomial, the p-bracket and the Potts model without 
changing the definition of the homology groups or chain spaces. This also raises the question of 
the relationship of the Khovanov homology with those constructions that have been made (e.g. 
G6l0 where the homology has been adjusted to fit directly with the dichromatic polynomial. We 
will take up this comparison in the next section. 

We now look more closely at the Potts model by writing a translation between the variables 
q, p and Q, v. We have 

-qp = Q~*v 

and 

q + q" 1 = QK 

and from this we conclude that 

q 2 - >/Qq + 1 = 0. 
Whence 

q = 

and 

1 = VQtVQ^ 

q~ 2 
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Thus 

„ -1 ±^1-4/0 

" = = " ( 2 ' 

For physical applications, Q is a positive integer greater than or equal to 2. Let us begin by 
analyzing the Potts model at criticality (see discussion above) where —pq = 1. Then 

1 -v/Q±v/Q^4 

p =~r — 2 — • 

For the Khovanov homology (its Euler characteristics) to appear directly in the partition function 
we want 

p=l. 

Thus we want 



Q± y/Q-A. 

Squaring both sides and collecting terms, we find that 4 — Q = tVQVQ ~ 4- Squaring once 
more, and collecting terms, we find that the only possibility for p = 1 is Q = 4. Returning to 
the equation for p, we see that this will be satisfied when we take a/4 = —2. This can be done 
in the parametrization, and then the partition function will have Khovanov topological terms. 
However, note that with this choice, q = —1 and so vj\[Q = —pq = 1 implies that v = —2. 
Thus e K — 1 = —2, and so 

e K = -l. 

From this we see that in order to have p = 1 at criticality, we need a 4- state Potts model with 
imaginary temperature variable K = (2n + l)z7r. It is worthwhile considering the Potts models at 
imaginary temperature values. For example the Lee- Yang Theorem ll20ll shows that under certain 
circumstances the zeros on the partition function are on the unit circle in the complex plane. 
We take the present calculation as an indication of the need for further investigation of the Potts 
model with real and complex values for its parameters. 

Now we go back and consider p = 1 without insisting on criticality. Then we have 1 = 
—vj (qy/Q) so that 

v = -qyQ = • 

From this we see that 



K _ _2-QtVQVQZA 

e i+v 2 

From this we get the following formulas for e A : For Q = 2 we have e K = ±i. For Q = 3, we 
have e K = ~ l± ^ 1 . For Q = 4 we have e A = — 1 . For Q > 4 it is easy to verify that e K is real and 
negative. Thus in all cases of p = 1 we find that the Potts model has complex temperature values. 
In a subsequent paper, we shall attempt to analyze the influence of the Khovanov homology at 
these complex values on the behaviour of the model for real temperatures. 
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4 The Potts Model and Stosic 's Categorification of the Dichro- 
matic Polynomial 

In ll26l Stosic gives a categorification for certain specializations of the dichromatic polynomial. 
In this section we describe this categorification, and discuss its relation with the Potts model. The 
reader should also note that the relations between dichromatic (and chromatic) homology and 
Khovanov homology are observed in [7] Theorem 24. In this paper we use the Stosic formulation 
for our analysis. 

For this purpose, we define, as in the previous section, the dichromatic polynomial through 
the formulas: 

Z[G](v, Q) = Z[G']{v, Q) + vZ[G"](v, Q) 
Z[*UG] = QZ[G}. 

where G' is the result of deleting an edge from G, while G" is the result of contracting that 
same edge so that its end-nodes have been collapsed to a single node. In the second equation, 
• represents a graph with one node and no edges, and • U G represents the disjoint union of 
the single-node graph with the graph G. The graph G is an arbitrary finite (multi-)graph. This 
formulation of the dichromatic polynomial reveals its origins as a generalization of the chromatic 
polynomial for a graph G. The case where v = — 1 is the chromatic polynonmial. In that case, 
the first equation asserts that the number of proper colorings of the nodes of G using Q colors is 
equal to the number of colorings of the deleted graph G' minus the number of colorings of the 
contracted graph G" . This statement is a tautology since a proper coloring demands that nodes 
connected by an edge are colored with distinct colors, whence the deleted graph allows all colors, 
while the contracted graph allows only colorings where the nodes at the original edge receive the 
same color. The difference is then equal to the number of colorings that are proper at the given 
edge. 

We reformulate this recursion for the dichromatic polynomial as follows: Instead of contract- 
ing an edge of the graph to a point in the second term of the formula, simply label that edge 
(say with the letter x) so that we know that it has been used in the recursion. For thinking of 
colorings from the set {1, 2, • • • , Q} when Q is a positive integer, regard an edge marked with x 
as indicating that the colors on its two nodes are the same. This rule conicides with our inter- 
pretation of the coloring polynomial in the last paragraph. Then G" in the deletion-contraction 
formula above denotes the labeling of the edge by the letter x. We then see that we can write the 
following formula for the dichormatic polynomial: 

Z[G] = J2 Q ]H] v e{H) 

HcG 

where H is a subgraph of G, \H\ is the number of components of H, and e(H) is the number of 
edges of H. The subgraphs H correspond to the graphs generated by the new interpretation of 
the deletion-contraction formula, where contraction is replaced by edge labelling. 
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Moving now in the direction of Euler characteristics, we let w = —v so that 

Z[G] = Z[G'] - wZ[G"} 
Z[m UG} = QZ[G] 

and 

Z[G] = J2 (-l) e{H) Q m w e{H) . 

HcG 

This suggests that differentials should increase the number of edges on the subgraphs, and that 
the terms Q^ H ^w e ^ should not change under the application of the (partial) differentials. Stosic's 
solution to this requirement is to take 

Q = q n 

and 

w = 1 + q + q 2 + ■ ■ • + q n 

so that 

Z[G] = (-l) e(/ V |H| (l + q + a 2 + • • • + q n ) e{H) . 
HcG 

To see how this works, we first rewrite this state sum (over states H that are sub-graphs of G) 
as a sum over enhanced states h where we define enhanced states h for a graph G to be labeled 
subgraphs h where a labeling of h consists in an assignment of one of the elements of the set 
S = {1, X, X 2 , ■ ■ ■ , X n } to each component of h. Regard the elements of S as generators of the 
ring R = Z[X]/(X n+1 ). Define the degree of X 1 by the formula deg(X l ) = n — i and let 

j(h) = n\h\ + deg(label(j)) 

iec(h) 

where the sum goes over all 7 in C(h), the set of components of h (each component is labeled 
from S and \h\ denotes the number of components in h). Then it is easy to see that 

Z[G] = J2 (-l) e{h) q j{h) 
heS(G) 

where S(G) denotes the set of enhanced states of G. 

We now define a chain complex for a corresponding homology theory. Let Ci(G) be the 
module generated by the enhanced states of G with i edges. Partial boundaries applied to an 
enhanced state h simply add new edges between nodes, or from a node to itself. If A and B arc 
components of an enhanced state that are joined by a partial boundary to form a new component 
C, then C is assigned the label when A and B have respective labels X % adn X° . This partial 
boundary does not change the labels on other components of h. It may happen that a component 
A is transformed by adding an edge to itself to form a new component A'. In this case, if A has 
label 1 we assign label X n to A' and otherwise take the partial boundary to be zero if the label of 
A is not equal to 1. It is then easy to check that the partial boundaries defined in this way preserve 
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j(h) as defined in the last paragraph, and are compatible so that the composition of boundary with 
itself is zero. We have described Stosic's homology theory for a specialization of the dichromatic 
polynomial. We have, as in the first section of this paper, 

Z[G] =Y^q J x(C' j (G)) =J2<l j x(H' j (G)) 

j j 

where x(C , - 7 (G)) denotes the complex defined above generated by enhanced states h with j = 
j(h), and correspondingly for the homology. 

Now lets turn to a discussion of the Stosic homology in relation to the Potts model. In the 
Potts model we have Q = q n is the number of spins in the model. Thus we can take any q 
so that q n is a natural number greater than or equal to 2. For example, we could take q to be 
an n-th root of 2 and then this would be a two-state Potts model. On the other hand, we have 
w = —v = 1 — e K as in our previous analysis for the Potts model. Thus we have 

-e K = q + q 2 H h q n . 

With q real and positive, we can take 

K = in + ln(q + q 2 + • • • + q n ), 

arriving at an imaginary temperature for the values of the Potts model where the partition func- 
tion is expressed in terms of the homology. If we take q — (1 + n) l l n , then the model will have 
Q = n + 1 states, and so, in this case, we can identify the enhanced states of this model as corre- 
sponding to the spin assignments of {1, X, X 2 , ■ ■ ■ , X n } to the subgraphs, interpreted as regions 
of constant spin. The partial boundaries for this homology theory describe particular (global) 
ways to transit between spin-labeled regions where the regions themselves change locally. Usu- 
ally, in thinking about the dynamics of a model in statistical physics one looks for evolutions that 
are strictly local. In the case of the partial differentials we change the configuration of regions at 
a single bond (edge in the graph) but we make a global change in the spin-labeling (for example, 
from X' 1 and X J on two separate regions to X J+J on the joined region). It is likely that the reason 
we see the results of this cohomology in the partition function only at imaginary temperature 
is related to this non-local structure. Nevertheless, the categorified homology is seen in direct 
relation to the Potts partition function and this connection deserves further examination. 



5 Imaginary Temperature, Real Time and Quantum Statistics 

The purpose of this section is to discuss the nature of imaginary temperature in the Potts model 
from the point of view of quantum mechanics. We have seen that for certain values of imaginary 
temperature, the Potts model can be expressed in terms of Euler characteristics of Khovanov 
homology. The suggests looking at the analytic continuation of the partition function to relate 
these complex values with real values of the temperature. However, it is also useful to consider 
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reformulating the models so that imaginary temperature is replaced with real time and the context 
of the models is shifted to quantum mechanics. To see how this works lets recall again the general 
form of a partition function in statistical mechanics. The partition function is given by the formula 

Z G (Q,T)=Y,e^ 1/kT)E ^ 

a 

where the sum runs over all states a of the physical system, T is the temperature, k is Boltzmann's 
constant and E(a) is the energy of the state a. In the Potts model the underlying structure of the 
physical system is modeled by a graph G, and the energy has the combinatorial form that we have 
discussed in previous sections. 

A quantum amplitude analogous to the partition function takes the form 

A G (Q,t) = Yl e(it/h)Eia) 

where t denotes the time parameter in the quantum model. We shall make precise the Hilbert 
space for this model below. But note that the correspondence of form between the amplitude 
Ag(Q, t) and the partition function Z G (Q, T) suggests that we make the substitution 

-1/kT = it/% 

or equivalently that 

t = (h/k)(l/iT). 

Time is, up to a factor of proportionality, inverse imaginary temperature. With this substitution, 
we see that when one evaluates the Potts model at imaginary temperature, it can be interpreted 
as an evaluation of a quantum amplitude at the corresponding time given by the formula above. 
Thus we obtain a quantum statistical interpretation of those places where the Potts model can be 
expressed directly in terms of Khovanov homology. In the process, we have given a quantum 
statistical interpretation of the Khovanov homology. 

To complete this section, we define the associated states and Hilbert space for the quantum 
amplitude A G (Q, t). Let H denote the vector space over the complex numbers with orthonormal 
basis {\cr)} where a runs over the states of the Potts model for the graph G. Define a unitary 
operator U (t) = e^*'^ 11 by the formula on the basis elements 

U{t)\a) = e {it/h)E{a) \a). 

The operator U (t) implicitly defines the Hamiltonian for this physical system. Let 

denote an initial state and note that 

U(t)\i/j) = Y^e {lt/h)E{a) \a) ) 



20 



and 

m(t)\i;) =Y / e Wh)E{ ° ) = A G (Q,t). 

Thus the Potts amplitude is the quantum mechanical amplitude for the state to evolve to the 
state U(t)\ip). With this we have given a quantum mechanical interpretation of the Potts model 
at imaginary temperature. 

Note that if, in the Potts model, we write v = e K — 1, then in the quantum model we would 
write 

Thus we can take t = —hKi, and if K is pure imaginary, then the time will be real in the quantum 
model. 

Now, returning to our results in Section 3 we recall that in the Potts model we have the 
following formulas for e K : For Q = 2 we have e K = ±i. For Q = 3, we have e K = ~ 1± ^ 1 . For 
Q = 4 we have e K = —1. It is at these values that we can interpret the Potts model in terms of 
a quantum model at a real time value. Thus we have these interpretations for Q = 2, t = hir/2; 
Q = 3, t = hn/6 and Q = 4, t — Hit. At these values the amplitude for the qunatum model is 
A G (Q, t) = J2a e {it/h)E{(T) and is given by the formula 

A G {Q,t)=Q N ' 2 {K{G)} 

where 

{K(G)} = Y,qi X (H-i(K(G))) 

j 

where H' 3 (K(G)) denotes the Khovanov homology of the link K(G) associated with the planar 
graph G and q — (1 — e lt / n ) / \[Q. At these special values the Potts partition function in its quan- 
tum form is expressed directly in terms of the Khovanov homology and is, up to normalization, 
an isotopy invariant of the link K(G). 



6 Quantum Statistics and the Jones Polynomial 

In this section we apply the point of view of the last section directly to the bracket polynomial. 
In keeping with the formalism of this paper we will use the bracket in the form 

<X> = (X> -<?<)<> 

with (O) = (q + OT 1 )- We have the formula for the bracket as a sum over enhanced states s : 

(K) = Y^(-l) nB{s) q j{s) 

s 
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where %(s) is the number of 5-type smoothings in s, A(s) is the number of loops in s labeled 1 
minus the number of loops labeled —1, and j(s) = tib(s) + A(s). In analogy to the last section, 
we define a Hilbert space TC (JC) with orthonormal basis { | s) } in 1 — 1 correspondence with the set 
of enhanced states of K. Then, for q = e l6 define the unitary transformation U : 7Y(/C) — > H(JC) 
by its action on the basis elements: 

U\s) = (-l) nB{s) q j ^\s). 
Setting \ip) — X) a l s )> we conclude that 

(k) = (mm- 

Thus we can express the value of the bracket polynonmial (and by normalization, the Jones 
polynomial) as a quantum amplitude when the polynomial variable is on the unit circle in the 
complex plane. 

There are a number of conclusions that we can draw from this formula. First of all, this 
formulation constitutes a quantum algorithm for the computation of the bracket polynomial (and 
hence the Jones polynomial) at any specialization where the variable is on the unit circle. We 
have defined a unitary transformation U and then shown that the bracket is an evaluation in the 
form (ip\U\ip). This evaluation can be computed via the Hadamard test II231 and this gives the 
desired quantum algorithm. Once the unitary transformation is given as a physical construction, 
the algorithm will be as efficient as any application of the Hadamard test. This algorithm requires 
an exponentially increasing complexity of construction for the associated unitary transformation, 
since the dimension of the Hilbert space is equal to the 2 c( - K ^ where c(K) is the number of cross- 
ings in the diagram K. Nevertheless, it is significant that the Jones polynomial can be formulated 
in such a direct way in terms of a quantum algorithm. By the same token, we can take the basic 
result of Khovanov homology that says that the bracket is a graded Euler characteristic of the 
Khovanov homology as telling us that we are taking a step in the direction of a quantum algo- 
rithm for the Khovanov homology itself. This will be the subject of a separate paper. For more 
information about quantum algorithms for the Jones polynonmial, see Ol [T4l US [25J . The form 
of this knot amplitude is also related to our research on quantum knots. See lfT6l . 
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